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Supersymmetric Construction of three-dimensional isospectral 

systems 

Yoshihide Yamada 

r-C ". Graduate School of Human and Environment Studies, 

Three-dimensional isospectral systems are constructed using the framework of supersym- 
metric quantum mechanics. In case the supercharge of first order in momentum is used, it is 
proved that the constructed systems reduce to a trivial generalization of one-dimensional sys- 
' terns. In case the second order supercharge is used, nontrivial three-dimensional isospectral 

systems are constructed. 

m . 

' j-h ' §1. Introduction 

+L> , The concept of super symmetry arose in the study of elementary particle physics. 

The term "supersymmetry" was originally used in reference to the symmetry of 
bosons and fermions, which behave according to different statistical laws. Super- 
symmetry is widely regarded as a necessary concept for unifying all the elementary 
forces. The concept of supersymmetry was originally formulated in the framework 
of field theory. However, it is not easy to analyze the results from supersymmetry 
(that is, whether or not a model exhibits the breaking of supersymmetry) in field the- 
ory. For this reason, supersymmetric quantum mechanics was proposed. 1 - 1 Although 
■ the original purpose of formulating supersymmetric quantum mechanics was to help 

obtain an understanding of supersymmetric quantum field theory, once this formu- 
lation of quantum mechanics was made, many interesting applications were found. 
Supersymmetric quantum mechanics is actually a reformulation of Darboux-Crum 
[ transformations and the factorization method. 2 ) Darboux-Crum transformations are 

techniques of generating a new Hamiltonian and its eigenfunctions from an existing 
Hamiltonian and its known eigenfunctions. The new spectrum is almost identical 
to the original. 3 ^' 4 ) The factorization method is a technique of solving Schrodinger 
equation by factorizing the Hamiltonian, which is the second order differential op- 
erator with respect to the spatial variables in position space, into the first order 
operators. 5 ' 

In quantum mechanics, supersymmetry refers to a relationship of isospectrality 
of two systems. More precisely, supersymmetric quantum mechanics is a formulation 
of quantum mechanics involving pairs of Hamiltonians (so called "partner Hamilto- 
nians") whose eigenvalues are identical (except for a possible ground state). The 
existence of such (non-trivial) pairs is possible because there is not a one-to-one cor- 
respondence between a spectrum and a potential. Indeed, for a given spectra there 
will be many potentials from which it can be derived. Supersymmetric quantum 
mechanics provides a method to construct pairs of isospectral systems. 6 ) 

An operator called intertwiner plays a central role in constructing isospectral 
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systems. Intertwiner was first introduced as factors of the Hamiltoinans, and there- 
fore they were of first order in the momentum. Later its generalization to the higher 
order in the momentum was proposed, which made it possible to construct other 
classes of isospectral systems. 7 )~ 9 ) 

However, most investigations of systems on super symmetric quantum mechan- 
ics and its higher order generalization are restricted to one-dimensional systems (or 
spherically symmetric three-dimensional systems). Two-dimensional generalization 
of supersymmetric quantum mechanics using first order intertwiner was proposed, 
which established an isospectrality of matrix Hamiltonian and a pair of scalar Hamil- 
tonians. 10 ) Two-dimensional isospectral systems without matrix Hamiltonian were 
constructed using the second order intertwiner. 11 )' 12 ) In this paper, we construct 
three-dimensional isospectral systems using supersymmetric quantum mechanics and 
its second order generalization. Two-dimensional isospectral systems based on first 
order intertwiner is also obtained from our result in three-dimensional case by taking 
the z-component of relevant quantity to zero. 

In section 2, we review the one-dimensional supersymmetric quantum mechanics 
and its second order generalization emphasizing the scheme for constructing isospec- 
tral systems. In section 3, we investigate three-dimensional generalization with first 
order intertwiner. In section4, we investigate three-dimensional generalization with 
second order intertwiner. 



§2. Review of one-dimensional case 



2.1. Supersymmetric quantum mechanics 

Let us see how one-dimensional isospectral systems are constructed using the 
framework of supersymmetric quantum mechanics. We consider a pair of nonrela- 
tivistic Hamiltonians; 

H = -d 2 + V(x), (2-1 
H = -d 2 + V(x), (2-2 
where d = d/dx. If these Hamiltonians satisfy the following relations; 

AH = HA, (2-3 
A*H = HA\ (2-4 

which are called "intertwining relations," then these systems are isospectral (except 
for zero mode of A or A^). This is because by multiplying A and A^ to 



respectively, we obtain 



Hi; E = E^ E , (2-5 

H$ E = (2-6 

H{A^ E ) = E(AiPe), (2-7 

H(A^ E ) = E{A^~). (2-8 
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From these relations, we see that H and H are basically isospectral (except that it 
is possible for the ground state of H or H to be missing), and their eigenfunctions 
are related as 

^ E oc Ai/> E , (2-9) 
ipE oc A^ E - (2-10) 

In the case the intertwiner of the form 

A = v{x)d + w(x) (2-11) 

is substituted to Eq. (12-3j) . we obtain the following equations for v,w,V and V; 

v' = 0, (2-12) 
v(V -V)- v" - 2w' = 0, (2-13) 
w(V — V) — vV' — w" = 0, (2-14) 

where prime denotes derivative with respect to the variable of the function. From 
Eq. (|2-12p we can set v = 1 without losing generality because A is determined only 
up to a multiplicative constant, as is evident in Eq. (|2-3p . Then from Eq. (|2-13p we 
obtain 

V-V = 2w', (2-15) 
and substituting this to Eq. (|2-14p . we find 

V = w 2 -w' + c, (2-16) 

V = w 2 + W > + C: (2-17) 

where c is an integration constant. The function w(x) remains arbitrary. By choosing 
the function w(x), isospectral systems V, V are obtained. The function w(x) is called 
the superpotential in supersymmetric quantum mechanics (though it corresponds 
to the derivative of the superpotential in supersymmetric quantum field theory). 
Products of A and A* are related to the Hamiltonians as follows, 

A^A = H-c, (2-18) 
AA^ = H-c. (2-19) 

As noted above, the ground state of H or H could be missing from the spectrum of 
the other. This is seen as follows: By defining the zero-mode functions of intertwiner 
A and A^ as 

A(jT = 0, (2-20) 
A ] ^ + = 0, (2-21) 
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we have 

H(f)~ = (A* A + c)(f>- = c(f>-, (2-22) 
H(p + = (AA ] + c)<p + = ccp + . (2-23) 

If (j)^ are regular and normalizable, they are eigenfunctions of H or H belonging to 
eigenvalue c which are ground states because from Eqs. (|2-18p and (|2-19p H — c and 
H — c are positive semi-definite^ — c > 0). Yet as (fx^ are written in terms of w(x) 
as 



exp( ± /^ W ), ,2-24) 



they are not normalizable simultaneously Thus if is normalizable, the 

ground state of H (H) is missing. And if neither of ^ are normalizable, spectra 
of H and H are identical and the missing of ground state does not occur (so called 
broken SUSY). The normalizability of ^ is decided from aymptotic behavior of 
superpotential w(x). 1 ^ 

By defining superhamiltonian H, supercharges Q and as 

»(V *•-«). 



we find the following supersymmetry algebra; 

[H,Q] = [H,Q*] = 0, (2-27) 

{Q,Q^} = H, (2-28) 

Q 2 = Qt2 = o, (2-29) 

where is the 2x2 null matrix. This is why this formulation is called supersymmetric 
quantum mechanics. 

2.2. Second order generalization 

Let us now consider a second order intertwiner; 13 )' 14 ) 

A = u{x)d 2 + v{x)d + w{x). (2-30) 
Substituting this to Eq. (|2-3|) . we obtain the following equations for u, v, w, V and V; 

v! = 0, (2-31) 

u(V -V)- u" - 2v' = 0, (2-32) 

v(V -V)- 2uV' - v" - 2w' = 0, (2-33) 

w(V -V)- uV" - vV - w" = 0. (2-34) 
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From Eq. (l2-31|) we can set u = 1 without losing generality for the same reason as 
the previous case. Then from Eq. (|2-32p we obtain 

V — V = 2v', (2-35) 

and substituting this to Eqs. (|2-33j) and (j2-34[) we have 

V + w = ^(v 2 - v') + c, (2-36) 
(V + w)" = 2(V + w)v' - 2v'V - vV, (2-37) 

where c is an integration constant. Eliminating (V + w) and multiplying v , we have 



Thus we obtain 



and from Eq. (|2-36|) 



-V^(^ + c, 2 -,V + ^-^Y. (2-38) 



v 2 , v" v' 2 d , nnn . 

V = — + c — v + - -j—r -\ — n > 2-39 

~ v 2 , v" v' 2 d , ndn . 

y = —r + c + v + - —r H — g , 2-40 

4 2u 4u 2 



v' v" v 12 d . A . 

4 2 2w Av 2 v 2 



where d is an integration constant. The function v{x) remains arbitrary. By choosing 
the function v(x) isospectral systems V,V are obtained. Products of A and are 
related to the Hamiltonians as follows, 

A* A = (H — c) 2 - d, (2-42) 
Ad = (H - c) 2 - d. (2-43) 

As there are two independent zero-mode functions of A, the spectra of H and H 
could be different at most two energy levels depending on their regularity and nor- 
malizability. This possible difference in spectra is rather delicate problem and it is 
usually analyzed based on Eqs. (|2-42p . (|2-43p and in terms of Darboux-Crum trans- 
formations. 8 )' 15 )' 16 ) 

The anticommutator of Q and becomes 

{Q, Qt} = H 2 - dl, (2-44) 



where 1 is the 2x2 identity matrix. This together with Eqs. (|2-27p and (|2 29|) form 
extended nonlinear superalgebra. 
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§3. Three-dimensional generalization with first order intertwiner 

Here we consider the three-dimensional generalization using the first order in- 
tertwiner. Hamiltonians are 

H = -A + V(r), (3-1) 
H = -A + V(r), (3-2) 

where r = (x, y, z), A = d^. + dy + d 2 z , and the intertwiner is the form 

A = v(r) ■ V + w{r), (3-3) 

where v{r) = (v x (r),v y (r),v z (r)). Substituting these to 

AH = HA, (3-4) 

we obtain the following equations for v,w,V and V; 

diVj + djVi = 0, (3-5) 
Av + 2Vw = (V - V)v, (3-6) 
v-VV + Aw = (V-V)w. (3-7) 

The general solution of Eq . (j3-5j) is 

v = a x r + b, (3-8) 

where a, b are arbitrary constant vectors. Thus the intertwiner has the form 

A = a- (r x V) + 6- V + w. (3-9) 

Substituting Eq. (|3-8jl into Eqs. (|3-6|) and (|3-7|) . these become 

2Vw = (V -V)(ax r + b), (3-10) 
a- (r x VV) + b- W + Aw = (V-V)w. (3-11) 

We find it convenient to deal with two cases, (i) a / and (ii) a = 0, separately. 

(i) When a / 0, we can take a to a unit vector without losing generality be- 
cause A is determined only up to a multiplicative constant. By making rotation of 
coordinates, we can set a = e z (the unit vector in the direction of the z axis); 

v = e z x (r - e z x b) + b z e z . (3-12) 

By translation of the origin r — e z x b — > r, this reduces to 

v = e z xr + b z e z , (3-13) 

which corresponds to a = e z and b = b z e z . Thus the intertwiner becomes 

A = e z -(rxV) + b z d z + w(r). (3-14) 
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Egs- tpnUP and (|34T|) become 

2Vw = (V-V)(e z x r + b z e z ), (3-15) 

e z ■ (r x W) + b z d z V + Aw = (V - V)w, (3-16) 

which reduces to the following equations using cylindrical coordinates (p,(p,z), 

2d p w = 0, (3-17) 

2d v w = (V- V)p 2 , (3-18) 

2d z w = (V-V)b z , (3-19) 



(S v + + ( ^ + -d p + -d% + d 2 z )w = (V- V)w. (3-20) 



When b z ^ 0, by differentiating Eqs. (|3-18p and (|3-19p with respect to p, we 
obtain a singular (on the z axix) potential 

V = V = V(p,0, t^b zH >-z. (3-21) 

For single- valuedness, V(p,£) must be a periodic function of £ with a period 27r6 z . 
Since H = H, the intertwiner becomes a symmetry operator of this system; 

[H,A)=0, A = d v + b z 8 z . (3-22) 

This is simply because [A, £] = 0. 

When 6^ = 0, Eq. (|3-19|) reduces to 9 2 io = 0, which together with Eq. (|3-17j) . 
yields 

w = w(ip). (3-23) 

Eq. (|3-18|) can be written as follows; 

V-V= 2 -^. (3-24) 
P 

By substituting Eas. iET^H and (ETMD to Ea. ([5^D|> . we have 

^ + ^M = ^M. (3-25) 

Thus we obtain 

^ = o l 3 ' 26 ) 

p 

where V pz is an arbitrary function of p, z, and 

A = d v + w((p). (3-28) 



8 



Y. Yamada 



And for single- valuedness, w(<p) must be a periodic function with a period 2ir. These 
potentials are singular on the z axis. To avoid the singularity, we have two means. 
One is to take w(ip) = const, which results in a trivial case V = V = V(p,z) and 
A = 8^ + const. The other is to consider that there is a cylindrical matter on the z- 
axis with the finite radius po , in such a case we may restrict the region of definition as 
p > Po an d assume a suitable boundary condition on wave function at p = po- In any 
way, Eqs, (|3-26p . (I3-27D and (I3-28D are trivial generalization of the one-dimensional 
case, because of the following reason: Firstly the following relation is satisfied by A 
in Eq.(GE28D; 

AH pz = H pz A, (3-29) 

where H pz = -p~ x d p pd p - d 2 z + V pz (p, z). Secondly, 

AH V = H V A, (3-30) 

where = —d 2 + w 2 ((p) — w'(ip) and Hp = —d 2 + w 2 (ip) + w'((p). By adding 
Eqs.flESl and d53^ devided by p 2 , we obtain 

A {y- + H » z ) = (^ + A - (3 ' 31) 

The potentials in H v /p 2 + H pz and H^/p 2 + H pz are of the form Eqs. (|3-26p and 



(ii) When a = 0, we can take b to a unit vector without losing generality. By 
making rotation of coordinates, we can set b = e x (the unit vector in the direction 
of the x axis). Then the intertwiner becomes 

A = d x + w(r), (3-32) 

and Eqs. ([3~T0j) and SSU\\ become 

2Vw = (V - V)e x , (3-33) 
d x V + Aw = (V - V)w. (3-34) 

Since y, z components of Eq. (|3-33p reduce to d y w = d z w = 0, we have 

w = w(x). (3-35) 

Then x component of Eq. (|3-33|) becomes 

V -V = 2w'(x). (3-36) 

By substituting Eqs. (|3-35p and (|3-36p to Eq. (|3-34p . we have 



d x V + w"(x) = 2w(x)w'{x). 



(3-37) 
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Thus we obtain 

V = w 2 (x)-w'(x) + V yz (y,z), (3-38) 

V = w 2 (x)+w'(x) + V yz (y,z), (3-39) 

where V yz (y,z) is an arbitrary function and 

A = d x + w(x). (3-40) 

This is a trivial generalization of one-dimensional case, because of the following 
reason: Firstly the following relation is satisfied by A in Eq. (|3-40p . 

AH yz = H yz A, (3-41) 

where H yz = —d 2 — d 2 + Vyz{y, z). Secondly, 

AH X = H X A, (3-42) 

where H x = —d 2 + w 2 {x) — w'(x) and H x = —d 2 + w 2 (x) + w'(x). By adding 
Eqs.(£RI]) and ([331]), we obtain 

A(H X + H yz ) = (H x + H yz )A. (3-43) 

By the results of (i) and (ii), we have proved that the three-dimensional gener- 
alization with the first order intertwiner reduces to a trivial generalization of one- 
dimensional case. 

§4. Three-dimensional generalization with second order intertwiner 

Here we consider the three-dimensional generalization using a second order in- 
tertwiner. Hamiltonians are Eqs. (|3-ip . ()3-2|) . and the intertwiner is the form 

A = giji^didj + v(r) ■ V + w(r). (4-1) 

Substituting these to Eq. (|3*4|) we obtain the following equations for gij,v,w,V and 
V; 

digjk + djg ki + d k g%j = 0, (4-2) 
Agij + diVj + d jVi = (V- V) 9ij , (4-3) 
2g ij d j V + Avi + 2d lW = (V- V)v u (4-4) 
g ij d i d j V + v-VV + Aw = (V-V)w. (4-5) 

The general solution of Eq. ()4-2p is 



5n =033^2 + a 22xj - 2a 23 x 2 x 3 - b u x 2 + 621x3 + cn, (4-6) 
522 =aiix§ + a 33 x\ - 2a 3 ix 3 a;i - &i 2 x 3 + b 32 xi + c 22 , (4-7) 
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933 =a.22xl + a n xl - 2a 12 x\X2 - h 3 xi + ^13^2 + c 33 , (4-8) 

512 = - 033^1^2 + «23^3^1 + ^31^2^3 ~ ^12^3 

, ^31 ^32 &11 - &22 / . n s 

+ ^x 1 -—x 2 X3 + C12, (4-9) 

523 = - CI11X2X3 + a 3 iX±X2 + 012X3X1 - <223X? 

1 bl2 ^13 ^22 1 (A tn\ 

+ -y x 2 - — x 3 - —X\ + c 23 , (4-10) 

531 = - 022^3^1 + ^12X2X3 + ^23X1X2 ~ 131^2 

1 ^23 hi .hi rA it\ 

+ -yx 3 - —X\ + —X 2 + C31, (4-11) 

where aij,bij,Cij are 20 arbitrary constant parameters. Then intertwiner Eq. ([4-l|) is 
rewritten as 

A = aijklj + bijkdj + Cijdidj + (<kjXi + Vj)dj + w, (4T2) 

by introducing (i) ij = ej^XjC^, (ii) 1221,013,032 so that ay satisfy Ojj = ay, (hi) 

633 so that S i=1 bij = 0, (iv) c 2 i,ci3,c 3 2 so that satisfy Cjj = and (v) dy = 

Ejfe=l a kk)5ij — CLij. 
@ 

We now consider a special solution of Eqs. (|4-3p . (|4-4j) and (j4- 5j) . We set a^j — 
bij = and Cij = diag(l, —1,0); 

A = Of - 8% + + w. (4-13) 

Then Eqs. (|4-3p . (|4-4p and (j4-5j) are rewritten as 

d i v j + d j v i = (V-V)c ij , (4-14) 

Cijdjiy + V) + 2d iW = (V- V)vi, (4-15) 

v- V{V + V) = 2(V-V)w. (4-16) 

By denoting x± = (xi ± X2)/V2, we can write the general solution of Eq. (|4-14p as 
follows; 

/3+0+) + /3_(x_) 

wi = x 3 yi H , (4-17) 

V2 = -X3V2 -j= , (4-18) 

v 3 = a — (cx + x^ + d + X- + d_x + ) (4-19) 

= a - (|x 2 + dixi) + {^x\ + d 2 x 2 ), (4-20) 

V -V = 2cx 3 + p' + (x+)+P'_(x.), (4-21) 

where a, c, d± = (di ±^2) A/2 are arbitrary constants, @±(x±) are arbitrary functions 
and 

2/1,2 = • (4-22) 
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By denoting y± = (y\ ± 2/ 2 )/v2, we can write the general solution of Eq. (|4-15p as 

~ 3 

V + V = 2V 3 {x 3 ) + -x\(y\ + y 2 _) + x 3 {y + P + + + a 2 (x 2 ) - ai(xi)) 
R 2 + B 2 

+ ^^- + 12{x2)- 71 (4-23) 
2w = c^x 2 . + + pL)x 3 + /?+/?- + 7i(^i) + 72(x 2 ). (4-24) 

In the above and 02(^2) must satisfy the condition 

aiixi) + a 2 (x 2 ) = v 3 (P' + + PL) - - (4-25) 

Now intertwiner becomes 



^ v^P'+ + PL)-y + P--y-P+ , 71+72 u 
+ g X3 + — 2 — ' ^ ' 

Here when v 3 = (.'. a = c = d-t = y± = 0), it reduces to a trivial generalization of 
two-dimensional case. Thus we consider the case v 3 7^ 0. Then Eq. (|4-16p is rewritten 
as four equations; 

2c{Acv 3 - 3(ca + d + d-)} = 2m 4 v 3 , (4-27) 
(4c + y+d+ + y-d-)a - (ca + d + d_)(P' + + p'_) = 2m 3 v 3 , (4-28) 
d + (p+a + y +7 ) + d_(/3_a + y_ 7 ) = « 3 {2m 2 + 3(y 2 + (4-29) 
d+(/3 +7 ) + 5_(/3_ 7 ) = t7 3 (2mi + + + a 2 - «i) (4-30) 

and 

Vs = + + ^x 2 + mix 3 , (4-31) 

where m are constants and a = a\(xi) + 02(2:2), 7 = 7i(xi) + 72(^2)- We found the 
general solution of equations for c 7^ as following; 

P± = h±y± + — , 
V± 

ai = -(V2htyl + a ), 

a 2 = V^hiyl + a , 
1 



71 = "4? 



„2, 2x .2 





(4-32) 




(4-33) 




(4-34) 




(4-35) 






£21 


(4-36) 


2/IJ ' 





72 = 4? { y 2 + (m2 " 4c2/l ?)^2 + 70 + 

where s, h± = (hi ± h 2 )/V2,q± = (qi db q 2 )j\[2 are arbitrary constants and 

2 hi(m 2 — Ac 2 h 2 ) 

7714 = 4c , m3 = 6v2c rti, mi = — ao, (4-37) 

v2 
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with conditions 

s±h,2 = s±q2 = s-qi = 0. (4-38) 

Transforming to the standard form as to X3 and omitting constants, we obtain the 
isospectral potentials; 

T/ 2 4, f 3 / 2 , 2\ , m 2 , 2.2] 2 , { ft 2 , 2 2 -> , 9l 1 

y = c X3 + + _ _3c h 1 jx 3 + ^^=(y + -y_) + — -c'jx 3 

m 2 3 l2 1 2 \ 2 , 2N , 1 /^ S l , S 2 



+ s? w + rf> + (jg - ?; + iM j + ri) + ^ $ + 3 

+ 2± (^^ + 2z k_±M^ + J =cAl , (4.39) 
4 4 y 2 ^ 

f> 2 4 , f 3 / 2 , 2N, m 2 o 2.2] 2 , / ^2 , 2 2 x , <?1 . 1 
1/4 4\ / /?71 2 3 j2 1 , o\ , 2 2n 1 f s l , s 2 

+ + vi) + - + -^l ) (y\ + </?) + ^ + g 

+ »-(« - *H + 1 
4 y2 4 y 2 ^ 

To repeat, c, m, h, q, s are constants which satisfy the conditions (|4-38p and y = cx+d. 
Thus we have constructed a class of new three-dimensional isospectral systems, which 
are forth order in X3. As these potentials have singularities of inverse square on 
= 0, we have to make restrictions on the region of definition as y± > 0. These 
systems have symmetry operator A^A and AA^ respectively by virtue of Eq. (|3-4|) 
and its Hermitian conjugate. Unlike one-dimensional case, A^A and AA^ are not 
written in terms of H and H in general. Thus the relationship between zero-modes 
of A and solutions of Schrodinger equation with Hamiltonian H is not clear. This 
makes the analysis of possible spectral difference more difficult. We leave it as a 
future problem. 



§5. Conclusion and Outlook 



Using the technique of supersymmetric quantum mechanics, we have constructed 
a class of three-dimensional isospectral systems. We have proved that the first or- 
der intertwiner case reduce to trivial generalizations of one-dimensional case. We 
have constructed new nontrivial three-dimensional systems using the second order 
intertwiner with special choice of "metric" gij = diag(l, — 1, 0). The constructed 
isospectral potentials do not contain arbitrary function of X3 but they are forth 
order in xj, (Of course as a trivial generalization of two-dimensional case, there 
are solutions which contain arbitrary function of X3). This limitation of functional 
form to the forth order happens also in the two-dimensional case with the met- 
ric §ij = diag(l, O). 17 ) We find it interesting that functions j3± are determined 
as Eq. (|4-32p because of nontrivial coupling to X3 component, though in the two- 
dimensional case with metric gij = (1,-1) there are wider classes of functions 
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(3±. These constructed potentials may not have direct physical significance for the 
present, but they might have some time. Whether these potentials may find direct 
physical application or not is left as a future problem. Even if there is no direct 
physical application, further analysis of these models will be meaningful because 
they can be used as a toy model for supersymmetric quantum field theory. For the 
same reason, the search for other isospectral systems based on different choices of 
metric gtj will be interesting. The problem of possible difference in spectra was left 
unsolved. In case of one-dimensional first order SUSY, the condition for missing of 
ground state is simple. Yet in case of one-dimensional second order SUSY, the sim- 
ple condition for spectral difference has not been found. Whether there is a simple 
condition or not is not sure. We might have to deal with it case by case depending 
on the models. In case of three-dimensional systems, the situation seems the same. 
It is an important problem to be studied for future works. In this paper, we have 
made a first step for application of supersymmetric quantum mechanics to construct 
three-dimensional isospectral systems having only one particle. The study of three- 
dimensional system is important because our physical space is three-dimensional, but 
in general, there are many particles concerning physical phenomenon. Whether the 
generalization of this method to the systems of more than two particles is possible 
or not is an open question. 
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